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Because most large galaxies contain a central black hole, and galaxies often merg#!, black- 
hole binaries are expected to be common in galactic nuclei. Although they cannot be imaged, 
periodicities in the light curves of quasars have been interpreted as evidence for binarieiP®, 
most recently in PG 1302-102, with a short rest-frame optical period of 4 yi®. If the orbital 
period matches this value, then for the range of estimated black hole masses the components 
would be separated by 0.007-0.017 pc, implying relativistic orbital speeds. There has been 
much debate over whether black hole orbits could be smaller than 1 ptP. Here we show that 
the amplitude and the sinusoid-like shape of the variability of PG 1302-102 can be fit by 
relativistic Doppler boosting of emission from a compact, steadily accreting, unequal-mass 
binary. We predict that brightness variations in the ultraviolet light curve track those in the 
optical, but with a 2-3 times larger amplitude. This prediction is relatively insensitive to the 
details of the emission process, and is consistent with archival UV data. Follow-up UV and 
optical observations in the next few years can test this prediction and confirm the existence 
of a binary black hole in the relativistic regime. 

Assuming PG 1302-102 is a binary, it is natural to attribute its optical emission to gas that 
is bound to each black hole, forming circumprimary and circumsecondary accretion flows. Such 
flows, forming “minidisks”, are generically found in high-resolution 2D and 3D hydrodynamical 
simulations that include the black holes in their simulated domaiiP®^. Assuming a circular orbit, 
the velocity of the lower-mass secondary black hole is 


/ 27r \ / GM \ 
\1 -f gj \4:7i‘^P J 


8,500 


1.5 


M 


l + qj VlO^-^Mg 


1/3 


P 


4.04 yr 


- 1/3 


km s 


-1 


or ~ 0.03c for the fiducial parameters above, where M = Mi -f M 2 is the total binary mass. Mi 2 
are the individual masses, q = M 2 /Mi < 1 is the mass ratio, P is the orbital period, and c is the 
speed of light. The primary’s orbital velocity is vi = qv 2 - Even if a minidisk has a steady intrinsic 
rest-frame luminosity, its apparent flux on Earth is modulated by relativistic Doppler beaming. The 
photon frequencies suffer relativistic Doppler shift by the factor D = [r(l — /9||)]“\ where T = 
(1 —is the Eorentz factor, /3 = u/c is the three-dimensional velocity v in units of the speed 
of light, and /3|| = /3 cos 0 sin i is the component of the velocity along the line of sight, with i and 0 
the orbital inclination and phase. Because the photon phase-space density oc Fu/y^ is invariant in 
special relativity, the apparent flux Fy at a fixed observed frequency v is modified from the flux of 
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a stationary source to "F°. The last step assumes an intrinsie power-law 

speetrum oc z/“. To first order mv/c, this eauses a sinusoidal modulation of the apparent flux 
along the orbit, by a fraetional amplitude AF^/Fi, = ±(3 — a){v cos(jy/c) sini Although light- 
travel time modulations appear at the same order, they are subdominant to the Doppler modulation. 
This modulation is analogous to periodie modulations from relativistie Doppler boost predietecP^ 
and observed for extrasolar planet^^H and for a double white dwarf binary!^, but here it has a 
mueh higher amplitude. 

The light-eurve of PG 1302-102 is well measured over approximately two periods (~ 10 
years). The amplitude of the variability is ±0.14 mag (measured in the optieal V bancP^I), eorre- 
sponding to AF^/Fi, = ±0.14. The speetrum of PG 1302-102 in and around the V band is well 
approximated by a double power-law, with a ^ 0.7 (between 0.50—0.55/xm) and a ^ 1.A (between 
0.55 — 0.6/im), apart from small deviations eaused by broad lines. We obtain an effeetive single 
slope ttopt = 1.1 over the entire V band. We eonelude that the 14% variability ean be attributed to 
relativistie beaming for a line-of-sight veloeity amplitude of v sin i = 0.074 c = 22, 000 km s“^. 

While large, this veloeity ean be realised for a massive (high M) but unequal-mass (low q) 
binary, whose orbit is viewed not too far from edge-on (high sin i). In Fig. 1, we show the required 
eombination of these three parameters that would produee a 0.14 mag variability in the sum-total of 
Doppler-shifted emission from the primary and the seeondary blaek hole. As the figure shows, the 
required mass is >10®'^ Mq, eonsistent with the high end of the range inferred for PG 1302-102. 
The orbital inelination ean be in the range of z = 60 — 90°. The mass ratio q has to be low q <0.3, 
whieh is eonsistent with expeetations based on eosmologieal galaxy merger modeld^, and also 
with the identifieation of the optieal and binary periods (for q >0.3, hydrodynamieal simulations 
would have predieted that the mass aeeretion rates fluetuate with a period several times longer than 
the orbital period^. 

As Fig. 1 shows, fully aeeounting for the observed optieal variability also requires that the 
bulk (/2 >80%) of the optieal emission arises from gas bound to the faster-moving seeondary blaek 
hole. We find that this eondition is naturally satisfied for unequal-mass blaek holes. Hydrodynam¬ 
ieal simulations have shown that in the mass ratio range 0.03 <q <0.1, the aeeretion rate onto the 
seeondary is a faetor of 10 — 20 higher than onto the primarjff^l. Beeause the seeondary eaptures 
most of the aeereting gas from the eireumbinary disk, the primary is “starved”, and radiates with 
a mueh lower effieieney. In the (M, q) ranges favoured by the beaming seenario, we find that the 
primary eontributes less than 1%, and the eireumbinary disk eontributes less than 20% to the to¬ 
tal luminosity, leaving the seeondary as the dominant souree of emission in the three-eomponent 
system (see details in Methods). 

The optieal light eurve of PG 1302-102 appears remarkably sinusoidal eompared to the best- 
studied previous quasi-periodie quasar binary blaek hole eandidate, whieh shows periodie burstd^. 
Nevertheless, the light eurve shape deviates from a pure sinusoid. In order to see if sueh devia¬ 
tions naturally arise within our model, we maximised the Bayesian likelihood over five parameters 
(period P, veloeity amplitude K, eeeentrieity e, argument of perieentre oj, and an arbitrary ref- 
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erence time to) of ^ Kepler orbiP^ and fit the observed optieal light-curve. In this procedure, we 
accounted for additional stochastic physical variability with a broken power-law power spectrum 
(i.e. a “damped random walk”^ described by two additional parameters. This analysis indeed re¬ 
turns a best-fit with a non-zero eccentricity of e = O.OQ^'^o og, although a Bayesian criterion does not 
favor this model over a pure sinusoid with fewer parameters (see Methods). We have considered an 
alternative model to explain PG 1302-102’s optical variability, in which the luminosity variations 
track the fluctuations in the mass accretion rate predicted in hydrodynamical simulation^^Hinillll]. 
However, the amplitude of these hydrodynamic fluctuations are large (order unity), and their shape 
is bursty, rather than sinusoid-likd ^^l^^l^^ , as a result, we find that they provide a poorer fit to the 
observations (see Fig. 2 and Methods). Furthermore, for mass ratios q >0.05, hydrodynamical 
simulations predict a characteristic pattern of periodicities at multiple frequencies, but an analysis 
of the periodogram of PG 1302-102 has not uncovered evidence for multiple peakP^. 

A simple observational test of relativistic beaming is possible due to the strong frequency- 
dependence of PG 1302-102’s spectral slope a = dlnF^/dlnu. PG 1302-102’s continuum spec¬ 
trum is nearly flat with a slope /3fuv = d\nFx/d\nX = 0 in the far UV (0.145-0.1525/xm) band, 
and shows a tilt with /Snuv = —0.95 in the 0.20-0.26/im near UV range (see Fig. 3 and Methods). 
These translate to apuv = —2 and cknuv = —1-05 in these bands, compared to the value a = 1.1 
in the optical. The UV emission can be attributed to the same minidisks responsible for the optical 
light, and would therefore share the same Doppler shifts in frequency. These Doppler shifts would 
translate into UV variability that is larger by a factor of (3 — a)Fuv/(3 — a)opt = 5/1.9 = 2.63 
and (3 — a)Nuv/(3 — a)opt = 4.05/1.9 = 2.13 compared to the optical, and reach the maximum 
amplitudes of ±37% (FUV) and 30% (NUV). 

PG 1302-102 has five separate UV spectra dated between 1992 and 2011, taken with instru¬ 
ments on the Hubble Space Telescope (HST) and on the GALEX satellite (see Fig. 3), as well as ad¬ 
ditional photometric observations with GALEX at 4 different times between 2006 and 2009 (shown 
in Fig. 2). The brightness variations in both the far- and near-UV bands show variability resem¬ 
bling the optical variability, but with a larger amplitude. Adopting the parameters of our best-fit 
sinusoid model, and allowing only the amplitude to vary, we find that the UV data yields the best- 
fit variability amplitudes of^Fy/Fy |fuv= ±(35.0 ±3.9)% and AFi./F,^ |nuv= ±(29.5 ±2.4)% 
(shown in Fig. 2). These amplitudes are factors of (2.57 ± 0.28) and (2.17 ± 0.17) higher than in 
the optical, in excellent agreement with the values 2.63 and 2.13 expected from the corresponding 
spectral slopes. 

Relativistic beaming provides a simple and robust explanation of PG 1302-102’s optical 
periodicity. The prediction that the larger UV variations should track the optical light-curve can be 
tested rigorously in the future with measurements of the optical and UV brightness at or near the 
same time, repeated two or more times, separated by a few months to ~ 2 years, covering up to 
half of the optical period. A positive result will constitute the first detection of relativistic massive 
black hole binary motion; it will also serve as a confirmation of the binary nature of PG 1302-102, 
remove the ambiguity in the orbital period, and tightly constrain the binary parameters to be close 
to those shown in Fig. 1. 
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Figure 1 | Binary parameters producing the optical flux variations of PG 1302-102 by 

relativistic boost. Shaded regions mark combinations of binary mass M, mass ratio 
q = M 2 /M 1 , and inclination i causing >13.5% flux variability (or line-of-sight velocity 
amplitude {v/c)smi > 0.07), computed from the Doppler factor with the effective 
spectral slope of a = l.l in the V band, including emission from the primary, as well as 
from the secondary black hole. The three panels assume fractions /2 = 1.0,0.95, or 0.8 
of the total luminosity arising from the secondary black hole; these values are consistent 
with fractions found in hydrodynamical simulation^^(seeA/er/zo(ij'). 
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MJD-49100 


Figure 2 | The optical and UV light-curve of PG 1302-102. Black points with la errors 

are optical dataP, superimposed with a best-fitting sinusoid (red dashed curve). The solid 
black curve shows the best-fit relativistic light-curve. The blue dashed curve shows the 
best-fit model obtained by scaling the mass-accretion rate found in a hydrodynamical 
simulation of an unequal-mass (g = 0.1) binarv*^. The additional circular data points with 
la errors show archival near-UV (red) and far-UV (blue) spectral observations; the red 
triangles show archival photometric near-UV data-points (see Fig. 3). The UV data include 
an arbitrary overall normalisation to match the mean optical brightness. The dotted red 
and blue curves show the best-fit relativistic optical light curve with amplitude scaled up 
by factors of 2.17 and 2.57, which best match the NUV and FUV data, respectively. 
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Figure 3 | Archival ultraviolet spectra of PG 1302-102 from 1992-2011. Far- and near- 

UV spectra obtained by the FOS and STIS instruments on the Flubble Space Telescope 
(FIST) and by GALEX are shown. Dates are in MJD (modified Julian day)-49100. Vertical 
yellow bands mark regions outside the spectroscopic range of both GALEX and FIST and 
contain no useful spectral data. From each spectrum, average flux measurements were 
computed in one or both of the two UV bands (shown in Fig. 2). GALEX photometric band 
shapes for FUV and NUV photometry are shown for reference as shaded blue and red 
curves, respectively. Additional GALEX NUV photometric data were also used in Fig. 2. 
The UV spectra show an offset by as as much as ±30%, close to the value expected from 
relativistic boost [see Methods). 
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METHODS 


V-band emission from a three-component system in PG 1302-102. Here we assume that the 
PG 1302-102 supermassive blaek hole (SMBH) binary system ineludes three distinct luminous 
components: a circumbinary disk (CBD), as well as an actively accreting primary and secondary 
SMBH. The optical brightness of each of the three components can be estimated once their accre¬ 
tion rates and the BH masses Mi and M 2 are specified. Using the absolute V-band magnitude of 
PG 1302-102, My = —25.81 and applying a bolometric correction BC ~ 10,^^ we infer the total 
bolometric luminosity of Lboi = 6.5(i?C'/10) x 10^® erg s“^. Bright quasars with the most massive 
SMBHs (M >10®M©), have a typical radiative efficiency of e = 0.3 1^. Adopting this value, the 

implied accretion rate is Mtot = -f^boi/(ec^) = 3 . 7 M 0 yr ^ 

We identify this as the total accretion rate through the CBD, and require that at small radii, 
the rate is split between the two black holes Mtot = M 2 -I- Mi with the ratio r] = M 2 /My 
Hydrodynamical simulation^ have found that the secondary captures the large majority of the 
gas, with 10 <rj <20 for 0.03 <q <0.1 (where q = M 2 /M 1 ). Defining the Eddington ratio of 
the disk as its accretion rate scaled by its Eddington-limited rate /i,Edd = ^i/^i,Edd with 
^Edd = -^Edd/O.lc^ (here -^Edd is the Eddington luminosity for the BH, and we have adopted 
the fiducial radiative efficiency of 0.1 to be consistent with the standard definition in the literature), 
we have 


/cBD,Edd « 0.068 (igtltle) • 


( 1 ) 

where the subscripts 1 and 2 refer to the primary and the secondary, and Mtot = Mi -f M 2 . We 
adopt a standard, radiatively efficient, geometrically thin, optically thick Shakura-Sunyaev (SS) 
disk modeP^ to compute the luminosities produced in the CBD and the circum-secondary disk 
(CSD). Although the secondary is accreting at a modestly super-Eddington rate, recent 3D radi¬ 
ation magneto-hydrodynamic simulations of super-Eddington accretion find radiative efficiencies 
comparable to the values in standard thin disk model^. On the other hand, the primary is accret¬ 
ing below the critical rate Mi < Madaf ~ 0. 027^/0. 3)^MEdd (with a the viscosity parameter) 
for which advection dominates the energy balance^. We therefore estimate its luminosity from a 
radiatively inefficient advection-dominated accretion flow (ADAE)|232ll, rather than a SS disk. This 
interpretation is supported by the fact that PG 1302-102 is known to be an extended radio source, 
with evidence for a jet and bends in the extended radio structureP^, features that are commonly 
associated with sub-Eddington sourced. 

Eor the radiatively efficient CBD and CSD, the frequency-dependent luminosity is given by 


/l,Edd — /cBD,Edd 


(1 + 9) 

l + Tj 


r j "IjEdd .. r, .7 

/2,Edd = V -- 1-37 

9 


- 0.0034 

/l,Edd j 

0.0034 ) 


( fcBB, 


Edd 


V 0.068 


1 + q 
1.05 


JL 

20 


0.05 


10 












integrating the loeal modified blaekbody flux over the area of the disk 
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where is the Planek funetion, is the frequeney-dependent absorption opaeity, and is 
the eleetron seattering opaeity. We eompute the radial disk photosphere temperature profile Tp by 
equating the viseous heating rate to the modified blaekbody flux 
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where a is the Stefan-Boltzmann eonstant and ks is Boltzmann’s eonstant. In solving for the 
photosphere temperature we work in the limit that k®® following Appendix A of ref. 

and we adopt rjg^o = QGM'^/c^ and {Ri^ Rout) = (2a, 200a), {rfsco^ for tho inner and 

outer radii of the CBD and CSD, respeetively. Here the superseript i refers to the disk, a is the 
binary separation, QGM is the loeation of the innermost stable eireular orbit for a Sehwarzshild 
blaek hole (our results are insensitive to this ehoise) and a(g/3)^/^ is the seeondary’s Hill radius 
(whieh provides an upper limit on the size of the CSEP^I). 


The optieal luminosity of an ADAF is sensitive to the assumed mierophysieal parameters 
and its eomputation is more eomplieated than for a thin disk. Here we first eompute a referenee 
thin-disk luminosity Lss for the primary, and multiply it by the ratio of the bolometrie luminosity 
of an ADAF to an equivalent thin-disk luminosity from ref. 


Ladaf 
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For ealeulating the referenee Lss, we adopted parameters eonsistent with ref. in partieular 
e = 0.1. Although the above ratio is for the bolometrie luminosities, we find that it agrees well 
with the faetor of 100 differenee in the V band shown in Figure 6 of ref. ^between ADAF and 
thin disk speetra with parameters similar to PG 1302-102 (IO^Mq, M = Madaf = 10“^'®MEdd, 
a ~ 0.3). 


Extended Data Fig. 1 shows the thin-disk CBD and CSD speetra for a total Eddington ratio of 
/cBD,Edd = 0.07, eonsistent with the high-mass estimates for PG 1302-102 needed for the beaming 
seenario (M = 10® “^ M© and q = 0.05). The red dot shows the redueed V-band luminosity of an 
ADAE onto the primary. The seeondary elearly dominates the total V band luminosity, with the 
primary eontributing less than 1%, and the CBD eontributing ~ 14%. In praetiee, the eontribution 
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from the CBD becomes non-negligible only for the smallest binary masses and lowest mass ratios 
(reaching 20% for M < IO^Mq and q < 0.025). 

We compute the contributions of each of the three components to the total luminosity, = 
L\ + L 2 + the corresponding total fractional modulation amplitude = 

(AL}^ + AL^)/-^tot’ for each value of the total mass M and mass ratio q. The primary is assumed 
to be Doppler-modulated with a line-of sight velocity vi = —qv 2 while the emission from the CBD 
is assumed constant over time (AL^g^ = 0). Extended Data Fig. 2 shows regions in the (M, q, i) 
parameter space where the total luminosity variation due to relativistic beaming exceeds 14%. This 
recreates Fig. 1 of the main text, but using the luminosity contributions computed self-consistently 
in the above model, rather than assuming a constant value of / 2 . Because the secondary is found to 
be dominant, the relativistic beaming scenario is consistent with a wide range of binary parameters. 


Model fitting to the PG 1302-102 optical light curve. We fit models to the observed light curve 
of PG 1302-102 by maximising the Bayesian likelihood C oc det|Cov^Cov^^|“^/^ exp(—x^/2), 
where 


X' = YT(Cov)-1Y, 


(5) 


and Y = O — M is the difference vector between the mean flux predicted in a model and the 
observed flux at each observation time U. Here Cov is the covariance matrix of flux uncertain¬ 
ties, allowing for correlations between fluxes measured at different fj. We include two types of 
uncertainties: (1) random (uncorrelated) measurement errors. 


Covg- = 


cr; 


i = j 
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( 6 ) 


where af is the variance in the photometric measurement for the data point (as reported in 
ref. 1^, and (2) correlated noise due to intrinsic quasar variability, with covariance between the 
and data points. 


Cov 


cr|,exp 


(1 -f z)td 


(7) 


The parameters ao and td determine the amplitude and rest-frame coherence time of correlated 
noise described by the damped random walk (DRW) modeP^, and the factor of (1 -f 2 ;) converts 
td to the observer’s frame where the are measured. These parameters enter the normalisation of 
the Bayesian likelihood, and this normalisation must therefore be included when maximising the 
likelihood over these parameter^. The covariance matrix for the total noise is given by Cov = 
Cov^ -f Cov^^. We assume both types of noise are Gaussian, which provides a good description 
of observed quasar variabilit)^^!. 


We then fit the following four different types of models to the data: 


• Relativistic beaming model with 5-i-2=7 model parameters: eccentricity, argument of peri- 
center, amplitude, phase and orbital period, plus the two noise parameters ao and To- 
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• Accretion rate model with 3+2=5 model parameters: amplitude, phase and period, plus the 
two noise parameters. This model assumes that PG 1302-102’s light-eurve traeks the mass 
aeeretion rates predieted in hydrodynamieal simulations. For near-equal mass binaries, sev¬ 
eral studies have found that the mass aeeretion rates fiuetuate periodieally, but they resemble 
a series of sharp bursts, unlike the smoother, sinusoid-like shape of PG 1302-102’s light- 
eurve. To our knowledge, only three studies to date have simulated unequal-mass (g < 0.1) 
SMBH binarie£3Iini]. The aeeretion rates for these binaries are less bursty; among all of 
the eases in these three studies, the q = 0.075 and q = 0.1 binaries in ref ^resemble 
PG 1302-102’s light-eurve most elosely (shown in Extended Data Fig. 3). Here we adopt 
the published aeeretion eurve for q = 0.1, and perform a fit to PG 1302-102 by allowing an 
arbitrary linear sealing in time and amplitude, as well as a shift in phase; this gives us the 
three free parameters for this model. (We find that the q = 0.075 ease provides a worse fit.) 

• Sinusoid model with 3+2 parameters: amplitude, phase and period, plus the two noise pa¬ 
rameters. This model is equivalent, to first order mv/c, to the beaming model restrieted to a 
eireular binary orbit. 

• Constant luminosity model with 2 parameters: for referenee only - to quantify how poor the 
fit is with only an amplitude plus the two noise parameters. 

In eaeh of these models above, we have fixed the mean flux to equal its value inferred from 
the optieal data. We have found that allowing the mean to be an additional free parameter did 
not ehange our results. The highest maximum likelihood is found for the beaming model, with 
best-fit values of (P = [IQQbjl'igg days, K = [0.065]1'!q qqq c, e = 0.09)')o;og, cosce = [—0.65]1 ')q gg, 
to = [718]^*134^ days), where the referenee point to is measured from MJD—49100. Uneertainties 
are eomputed with the emcee eodeP^, whieh implements a Markov Chain Monte Carlo algorithm, 
and whieh we use to sample the 7D posterior probability of the model given the G15 data. In 
praetiee, we employ 28 individual ehains to sample the posterior for 1024 steps eaeh. Throwing 
away the first 600 steps (‘burning in’), we run for 424 steps and for eaeh parameter, we quote best- 
fit values eorresponding to the maximum posterior probability, with errors given by the 85**^ and 
15**^ pereentile values (marginalized over the other six parameters) . The best-fit noise parameters 
are found to be (ud, td) = (0.049^° gg® mag, 37.6^1^5^ days). The best-fit model has a redueed 
X^/{N — 1 — 7) ~ 2.1, where N = 245 is the number of data points. 

To assess whieh of the above model is favoured by the data, we use the Bayesian Information 
Criterion (BIG), a standard method for eomparing different models, penalising models with a larger 
number of free parameter^. Speeifieally, BIC= —2\n.C + k\n.N, where the first term is evaluated 
using the best-fit parameters in eaeh of the models and where k is the number of model parameters. 
We find the following differenees ABIC between pairs of models: 

BICacc ~ BICBea^m = 4.0 (Beaming model preferred over accretion model) 

BICacc ~ BlCsin = 14.9 (Sin model strongly preferred over accretion model) 

BICsin — BICBesLia = —10.9 (Sin model strongly preferred over eccentric beaming model) 
BICconst — BICBeam = 11-5 (Beaming model strongly preferred over pure noise) 

BICconst — BICsin = 22.4 (Sin model strongly preferred over pure noise). 
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We conclude that a sinusoid, or equivalently the beaming model restriced to a circular binary, 
is the preferred model. In particular, this model is very strongly favoured over the best-fit accretion 
models (see Extended Data Fig. 3), with ABIC > 14.9. For the assumed Gaussian distributions, 
this corresponds to an approximate likelihood ratio of exp(—14.9/2) 5.7 x 10“^. Although 

our best-fit beaming model has a small non-zero eccentricity, the 7-parameter eccentric model is 
disfavoured (by ABIC = 11.5) over the 5-parameter circular case. 

We have conservatively allowed the amplitude of accretion rate fluctuations to be a free 
parameter in the accretion models, but we note that the accretion rate variability measured in 
hydrodynamic simulations exhibits large (order unity) deviations from the mean, even for 0.05 < 
g < 0.1 binaries EEHHl. in the accretion rate models, an additional physical mechanism needs to 
be invoked to damp the fluctuations to the smaller ~ 14% amplitude seen in PG 1302-102 (such as 
a more significant contribution from the CBD and/or the primary). 

Disk Precession. The lowest BIG model, with a steady accretion rate and a relativistic boost from 
a circular orbit, has a reduced = 2.1, indicating that the relativistic boost model with intrin¬ 
sic noise does not fully describe the observed light-curve. The residuals could be explained by 
a lower-amplitude periodic modulation in the mass accretion rate, which is expected to have a 
non-sinusoidal shape (i.e. with sharper peaks and broader troughs, as mentioned abovef^. Alter¬ 
natively, the minidisks, which we have implicitly assumed to be co-planar with the binary orbit, 
could instead have a significant tilP^. 

A circum-secondary minidisk that is tilted with respect to the binary’s orbital plane will 
precess around the binary angular momentum vector, causing additional photometric variations due 
to the changing projected area of the disk on the sky. The precession timescale can be estimated 
from the total angular momentum of the secondary disk and the torque exerted on it by the primary 
black hole. The ratio of the precession period to the binary’s orbital period is 

-Pprec _ 8 -sjl -\- q 

-Porb a/s cos (5 

where we have chosen the outer edge of the minidisk to coincide with the secondary’s Hill sphere 
= (g/3)^Pa, for binary semi-major axis a. This choice gives the largest secondary disk and the 
shortest precession rates. The angle 5 between the disk angular momentum vector and the binary 
angular momentum vector can range from —7r/2 to 71/2. For small binary mass ratios, consistent 
with the relativistic beaming scenario, the precession can be as short as 4.8Porb, causing variations 
on a timescale spanning the current observations of PG 1302-102. The precession timescale would 
be longer (> 20Porb) for a smaller secondary disk tidally truncated at 0.27g°-^a® and with a more 
inclined (45°) disk. 

Archival UV data. FUV (0.14-0.175/rm) and NUV (0.19-0.27/im) spectra of PG 1302-102 were 
obtained by the Hubble Space Telescope (HST) and the Galaxy Evolution Explorer (GAFEX) 
since 1992. HST/Faint Object Spectrograph (FOS) NUV spectra were obtained on July 17, 1992 
(pre-COSTARj®. HST/Space Telescope Imaging Spectrograph (STIS) FUV spectra were obtained 
on August 21, 2001^. GAFEX FUV and NUV spectra were obtained on March 8, 2008 and April 
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6, 2009 and HST/Cosmic Origins Spectrograph (COS) FUV spectra were obtained on January 28, 
2011. All data are publicly available through the Mikulski Archive for Space Telescopes (MAST) 
at archive.stsci.edu. All measurements have been spectrophotometrically calibrated, and binned or 
smoothed to l-SA resolution. The spectra (shown in Fig. 2 in the main text) have errors per bin 
typically less than 2% and published absolute photometric accuracies are better than 5%. 

From each spectrum, average flux measurements (shown in Fig. 2 in the main text) were 
obtained in one or both of two discrete bands: FUV continuum from 0.145-0.1525/im for FUV 
(a range chosen to avoid the Lya line) and NUV continuum from 0.20-0.26/im. For the GALEX 
NUV photometric data (also used in Fig. 2) we adopted a small correction (0.005 mag) for the 
transformation from the GALEX NUV to our NUV continuum band. GALEX EUV photometric 
data were not used because of the significant contribution from redshifted Lya. Note that the broad 
lines in the UV spectra (in Eig. 3) do not show a large AA = {v/c)\ ~ 140A Doppler shift. This is 
unsurprising, since the broad line widths (2,500-4,500 km s“^) are much smaller than the inferred 
relativistic line-of-sight velocities, and are expected to be produced by gas at larger radii, unrelated 
to the rapidly orbiting minidisks producing the featureless thermal continuum emissiorP^. 
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Extended Data Figure 1 | Model spectrum of PG 1302-102. Circumbinary (dashed blue) and 
circumsecondary (solid black) disk spectra for a total binary mass of 10®'^, binary mass ratio q = 
0.05, and ratio of accretion rates M 2 /Mi = 20. A vertical dashed line marks the center of the 
V-band and the approximate flux from an advection-dominated accretion flow (ADAF) is shown 
as a red dot for the V-band contribution of the primary. The spectrum for a radiatively efficient, 
thin disk around the primary is shown by the thin red dashed curve for reference. 
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i (From Face-On) 



Extended Data Figure 2 | Parameter combinations for which the combined V-band luminos¬ 
ity of the three-component system varies by the required 0.14 mag. M is the binary mass, q 
is the mass ratio, and i is the orbital inelination angle. This figure is analogous to Fig. 1, exeept 
instead of adopting a fraetional luminosity eontribution /2 by the seeondary, the luminosities of 
eaeh of the three eomponents are eomputed from a model: the primary’s luminosity is assumed to 
arise from an ADAF, while the seeondary’s luminosity is generated by a modestly super-Eddington 
thin disk. Emission from the eireumbinary disk is also from a thin disk, and is negligible exeept 
for binaries with the lowest mass ratio g <0.01 (see text). 
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Extended Data Figure 3 | Model fits to PG 1302-102’s optical light curve. Best-fitting curves 
assuming relativistic boost from a circular binary (solid black curves), a pure sinusoid (red dotted 
curves) and accretion rate variability adopted from hydrodynamical simulation^ (blue dashed 
curves) for a g = 0.075 (a) and a g = 0.1 (b) mass-ratio binary. The grey points with la errors 
bars show the data for PG 1302- 102P. 


19 



















